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We verified operational approach based on direct measurement of entanglement degree for bipar- 
tite systems. In particular spectral distributions of single counts and coincidence for pure biphoton 
states generated by train of short pump pulses have been measured and entanglement quantifier 
calculated. The approach gives upper bound of entanglement stored in total biphoton states, which 
can reach extremely high value up to 10 4 — 10 5 . 
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Recently entanglement became from an argument of 
debate on foundations on quantum mechanics a re- 
source for developing quantum technologies (as quan- 
tum communication, q-calculus, q-imaging, q-metrology, 
etc.). In this sense a precise and easy-implementable 
method for characterizing entanglement properties is a 
fundamental tool. There are several measures that quan- 
tify entanglement of quantum bipartite states both in 
discrete and continues variables [2j. The most popular 
ones are Schmidt rank, entropy, concurrence, etc. From 
a theoretical point of view all these quantifiers can be 
easily evaluated from the density matrix describing the 
given state, thus there is no problem to find the meaning 
of a chosen entanglement measure for a certain density 
matrix. From an experimental side and following the 
same paradigm it would be necessary to perform quan- 
tum tomography procedure (either complete or reduced) 
over the state [3| under study and then to extract the 
entanglement measure from reconstructed density ma- 
trix. However this scheme is not optimal in the case 
of high-dimensional systems when number of measure- 
ments for complete reconstruction of the density matrix 
grows quadratically with dimension of Hilbcrt space. An 
alternative is exploiting particular links between entan- 
glement and measurable characteristics of the state, just 
measuring some auxiliary parameters of the state. In par- 
ticular the parameter R q (Fedorov's ratio Q), defined as 
the ratio of the single-particle and coincidence distribu- 
tions widths in q-space, can be rather easily measured 
in contrast to all other entanglement quantifiers. R can 
be understood qualitatively from entropy approach to en- 
tanglement: high entanglement leads to the better knowl- 
edge about composite bipartite system (narrower coinci- 
dence distribution) and the worse knowledge about indi- 
vidual subsystcm(s) (wider single-particle distribution). 



Both clear physical meaning and opcrationability return 
R to an extremely useful tool for entanglement control. 
It has been proved 0] that for double-Gaussian bipar- 
tite states the parameter R coincides exactly with the 
Schmidt number K . Moreover their values remain quite 
close even for special classes of non-double-Gaussian wave 
functions like those describing Spontaneous Parametric 
Down-Conversion (SPDC). This approach can be used 
both for continues and discrete variables, so the param- 
eter R has rather universal sense: earlier it has been 
successfully applied for exposure strong entanglement 
anisotropy in spatial distributions of biphotons [5[ . 

The present paper is devoted to verification of the oper- 
ational entanglement quantifier R for two-photon states 
entangled in frequency domain as an efficient alternative 
to other quantifiers (for example in specific experiments 
visibility of interference pattern can provide some knowl- 
edge about entanglement Q). Indeed such states be- 
long to multi-dimensional Hilbcrt space and can posses 
extremely high entanglement degree (up to several hun- 
dreds) that makes its very perspective objects for quan- 
tum information and quantum communication. 

Let us remind first the main features of the biphoton 
spectral wave function for SPDC with the type-I degen- 
erate collinear phase matching and with a pulsed pump. 
The main general expression has the form 3, S] 
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where r is the pump-pulse duration, L is the length of 
the crystal, vi and v 2 are deviations of frequencies of the 
signal and idler photons ui 2 from the central frequencies 
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Up/ 2, 2 1 "C Up, u p is the central fre- 
quency of the pump spectrum, A and B are the temporal 
walk-off and dispersion constants 
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are the group velocities of the pump and 
ordinary waves, and k\ and k p are the wave vectors of 
signal and pump photons. 

The wave function {T]) can be served for determining 
coincidence and single-particle biphoton spectra. These 
spectra are significantly different in the cases of short an 
long pump pulses. The control parameter separating the 
regions of short and long pulses is given by 
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i.e., it is equal to the ratio between the double pump- 
pulse duration and the difference of times required for 
the the pump and idler/signal photons for traversing all 
the crystal. Pump pulses are short if rj <C 1 and long if 
t] 3> 1 and, typically, r\ ~ 1 at r ~ lps. 

In the most interesting case of short pump pulses the 
FWHM of the coincidence and single-particle spectra 
were found analytically to be given by 
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The degree of entanglement of the state (fTJ) was charac- 
terized by two parameters, the Schmidt number K and 
the parameter R defined as the ratio of the single- to 
coincidence spectral widths of the corresponding photon 
distributions, R^ = Au s /Au c AA S /AA C = R\ Q. 
Both parameters were calculated analytically as function 




FIG. 1: Quantifiers calculated analytically (R^ Q) and mirni 
ically (K for LiI03 crystal and the pump wavelength \ p 
400 nm. 



close to each other. Also, for the same wave function 
(P), K{t) was calculated numerically and the result 
was found to be in a very good agreement with the an- 
alytical function R(t) of Ref. Q, Fig. [Q The de gree of 
spectral entanglement reaches rather high values in the 
whole range of pump-pulse durations and especially high 
in the case of sufficiently short and long pulses (far from 
the minimum that localizes at 77 ~ I or r ~ 1 ps) . 

Evidently the parameter ([3|) accumulates factors to be 
rather easily controlled in experiment, namely length of 
the crystal and pump pulse duration. Also the longitu- 
dinal walk-off effect can be evolved via A by choosing a 
crystal with appropriate dispersion properties as a source 
of biphotons. 

The experiment(Fig.2) was performed in fs pulsed 
regime using a Mode-Locked Titanium-Sapphire laser at 
a working wavelength of Xm = (795.0 ± 0.1) nm with a 
AA/_r = (5.9 ± 0.1) nm. After doubling in frequency to 
Xp = (397.5 ± 0.2) nm and AA P = (1.8 ± 0.1) nm, cor- 
responding to a pulse duration of r = (186 ± 30) fs, the 
pump beam is addressed to a L1IO3 crystal where type 
I collinear SPDC is produced. After eliminating the UV 
pump, biphotons are split on a beam-splitter and fed to 
two photodetection apparatuses (consisting of red glass 
filters and SPAD detectors). In front of each detector it 
is placed a monochromator with a variable spectral reso- 
lution. Light is focused by a lens (/=20 cm) inside each 
monochromator and a couple of lenses with the same 
focal distance are placed on the output of the monochro- 
mators for focalizing the output packet spread by all the 
optics of the apparatus and by the slits. 
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FIG. 2: Sketch of the experimental set up: a duplicated ti- 
tanium sapphire laser beam pumps a L1IO3 crystal produc- 
ing collinear SPDC. After a Half Wave Plate (HWP) the 
biphotons are split on a beam-splitter (BS) and fed to SPAD 
detectors, whose output feeds Time-to Amplitude Converter 
(TAC), counters and Multi Channel Analyzer (MCA). 



of the pump-pulse duration r and were found to be very 



Following to general idea of the i?-quantifier method 
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all experiments were performed as a consecution of cor- 
responding measurements of spectral distributions both 
in single counts and coincidences. To measure the co- 
incidence distribution one of the two monochromators 
(by convention selecting idler wavelength) is fixed at the 
central wavelength of SPDC (795 nm) whilst the other 
(signal) scans in a range around this value. In order to 
check R dependence on the control parameter 77 the mea- 
surements have been repeated for crystals with different 
length (but with the same orientation), namely L=10 
mm and L~5 mm. 

First of all we have performed measurements with a 10 
mm L1IO3 crystal. To measure the coincidence distribu- 
tion we have used monochromators with 0.2 nm resolu- 
tion. A narrow peak with AA C = (0.29 ± 0.03) nm has 
been obtained (Fig. 3) while scanning the wavelength in 
signal channel. It is clearly seen that the pump width is 
6.2 times larger than the coincidence spectrum. 
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FIG. 3: Normalized coincidence distribution for 10 mm length 
L1IO3 sample. Points with error bars correspond to the measure- 
ments while the solid line shows the fitting by a Gaussian function 
with FWHM =(0.29 ± 0.03)nm. 

The spectral distribution of single counts was per- 
formed using a monochormator in the transmission arm 
with a spectral resolution of 1 nm, finding a single counts 
spectral width AA S = (101±1) nm (Fig. 4). The asymme- 
try of the right wing in the measured spectrum might be 
caused by falling spectral sensitivity of monochromator 
in the long-wavelength range. 

Thus the distribution of single counts is 56 times 
larger than that of the pump, corresponding to an ex- 
perimental ratio between widths of the two distributions 
R u = (349 ± 43). The complete results are collected in 
the table I. 

Then the same measurement has been performed with 
a crystal of 5 mm length obtaining AA S = (115 ± 1) nm 
and AA C = (0.64 ± 0.06) nm, corresponding to R u = 
(179 ± 18). Considering that the pulse duration of the 
pump did not change we can observe that in this case the 
pump is 2.8 times greater than the coincidence spectrum 
whereas the single counts distribution is 64 times larger 
than the width of the pump. Corresponding results are 
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FIG. 4: Normalized single counts distribution for 10 mm length 
L1IO3 sample. 

TABLE I: Results for 10 mm LHO3 sample. The pump pa- 
rameters t and AA P are measured independently by using a 
cross-correlator (with first harmonic) and a spectrometer. 



Theory Model 


Experimental results 


Pump 


A p = 397. 5nm 


A p = 397.5 ± 0. 2nm 


AA P = 1.2571m 


AAp = (1.8 ±0.1)71771 


r = 186/s 


t — 186 i 30/s 


Coincidence Distribution 


A c = 79571771 


A c = 795.0 ± 0.271771 


AA C = 0.3271771 


AA C = (0.29 ± 0.03) 


Single Counts Distribution 


AA S = 10071771 


AA S = (101 ± l)nm 


i?-Quantifier 


Rui-Rx- 316 


R^~ Rx~ (349 ± 43) 



listed in the table II. 



TABLE II: Results for 5 mm LiIQ 3 sample 



Theory Model 


Experimental results 


Coincidence Distribution 


AA C = 0.6371771 


AA C = (0.64 ± 0.06) 


Single Counts Distribution 


AA S = 10071771 


AA S = (115 ± l)nm 


i?-Quantifier 


Rui ~ R\ ~ 158 


R^~ Rx~ (179 ± 18) 



Certainly the most comprehensive contribution to the 
experimental verification of the i?-quantificr validity in 
frequency domain would be testing dependence R on 
the control parameter rj through all available parameters. 
Agreement between the measured and predicted shapes 
of the curve shown at Fig.l would confirm completely the 
adequacy of the i?-quantifier approach. Looking at (|3j [4j) 
it is seen that the " working" parameters might be r and 
L. However it is difficult to change the pump pulse dura- 
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tion in a wide range with the same laser (from hundreds 
of fs to dozens of ps) and pass through the minimum at 
the curve R{ji) varying this parameter only. That is why 
we took into account the fact of linear dependence of R 
on the sample length L through the coincidence distri- 
bution width (UJ). So doubling the sample length leads 
to doubling the entanglement degree for the short pump 
pulse regime of SPDC. This fact is clearly illustrated by 
the obtained results, demonstrating the validity of the 
conceptual scheme even if, unfortunately, subsequent de- 
crease of the sample length reduces the biphoton flux too 
much for measuring widths AA S and reaching the mini- 
mum of R(rf). 

One should mention that in some works (see, e.g., 
[lot ) the 2D- photon distribution was measured and plot- 
ted in the (wi,^) plane. Surprisedly that such 2D- 
distributions are very useful for making entanglement 
analysis of biphoton states in context of conditional pure 
state preparation, dispersion broadening [Til . [l3| . etc. 
This would allow evaluating i?-quantificr directly. We 
would like also to mention that a measurement of spectral 
width could also be obtained by transforming spectrum 
in time difference through a fiber 11 1 . 

Finally the question arises "how much entanglement 
can be stored in a given state". We remind that the 
method of measurement of R in frequency domain im- 
plies that other degrees of freedom of biphotons should be 
fixed. First off all we mean spatial distribution of photon 
pairs emitted by crystal, which must be illuminated by 
selecting particular narrow-angular segments with small 
pinholes. Vice versa, measurement of i?-quantifier in spa- 
tial variables (Rangie) implies fixing frequencies of both 
photons, which is usually performed with narrow-band 
filters |5j. Therefore the total entanglement potentially 
stored in the SPDC state can be conventionally accepted 
as a product R an gie * Ru ■ On the one hand this is a good 
convention since quantifier R has the physical meaning of 
effective Hilbert space, so the total dimension of biphoton 
states in joint polarization-frequency-spatial space would 



be R t 



2 * Rangie * Ru ■ However we suggest to use 



this just as an upper bound of the stored entanglement 
taking into account the existing links between frequency 
and spatial degrees of freedom caused by phase match- 
ing. Speaking in terms of numbers this value can reach 
really huge amount R tot w 10 4 - 10 5 [3]! 

In conclusion we have verified the approach of Ref. 
[H addressed to estimate the entanglement degree in 



frequency domain for pure biphoton states generated 
by short pump pulses. An approach that can find 
widespread application in quantum technologies exploit- 
ing entanglement of biphotons. The generalization to the 
mixed two-photon states will be discussed elsewhere. 
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